Valence-shell photoionization of chlorine-like Ar+ ions 

A. M. Covington, A. AguilarQi. R. Covington, G. Hinojosa0C. A. Shirley, and R. A. Phaneuf 
Department of Physics, MS 220, University of Nevada, Reno, NV 89557-0058, USA 

I. Alvarez and C. Cisneros 

Instituto de Ciencias Fisicas, Universidad Nacional Autonoma de Mexico, 

Apartado Postal 6-96, Cuernavaca 62131, Morelos, Mexico. 

I. Dominguez-Lopezj^ M. M. Sant'Anna|§] and A. S. Schlachter 
Advanced Light Source, Lawrence Berkeley National Laboratory, 1 Cyclotron Road, Berkeley, CA 94720, USA 



o 



O 
(N 

m 

6 

u 



C. p. Ballance 
Department of Physics, Auburn University, Auburn, AL 36840, USA 

B. M. McLaughlii]^ 

Institute for Theoretical Atomic and Molecular Physics, Harvard Smithsonian Center for Astrophysics, 

60 Garden Street, MS- 14, Cambridge, MA 02138, USA 

(Dated: February 27, 2013) 

Absolute cross-section measurements for valence-shell photoionization of Ar"'" ions are reported 
for photon energies ranging from 27.4 eV to 60.0 eV. The data, taken by merging beams of ions 
and synchrotron radiation at a photon energy resolution of 10 meV, indicate that the primary ion 
beam was a statistically weighted mixture of the ^P^/2 ground state and the ^Pi/2 metastable 
state of Ar+. Photoionization of this C^-like ion is characterized by multiple Rydberg series of 
autoionizing resonances superimposed on a direct photoionization continuum. Observed resonance 
lineshapes indicate interference between indirect and direct photoionization channels. Resonance 
features are spectroscopically assigned and their energies and quantum defects are tabulated. The 
measurements are satisfactorily reproduced by theoretical calculations based on an intermediate 
coupling semi-relativistic Breit-Pauli approximation. 
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I. INTRODUCTION 



Photoionization (PI) of ions is a fundamental process 
of importance in many high-temperature plasma envi- 
ronments, such as those occurring in stars and nebu- 
lae [T and in inertial-confinement fusion experiments [2J . 
Quantitative measurements of photoionization cross sec- 
tions for ions provide precision data on ionic structure 
and guidance to the development of theoretical models 
of multi-electron interactions. In addition, the opacity 



* Present address: Advanced Light Source, Lawrence Berkeley Na- 
tional Laboratory, 1 Cyclotron Road, MS 7-100, Berkeley, CA 
94720, USA 

tPresent address: Instituto de Ciencias Fisicas, Universidad Na- 
cional Autonoma de Mexico. Apartado Postal 6-96, Cuernavaca 
62131, Morelos, Mexico. 

-f Present address: Centre Nacional de Metrologia, Queretaro, Qro. 
76900, Mexico 

S Present address: Instituto de Fisica, Universidade Federal do Rio 
de Janeiro, Caixa Postal 68528, 21945-970 Rio de Janeiro RJ, Brazil 
^Present address: Centre for Theoretical Atomic, Molecular and 
Optical Physics (CTAMOP), School of Mathematics and Physics, 
The David Bates Building, 7 College Park, Queen's University of 
Belfast, Belfast BT7 INN, United Kingdom 



databases [3H5] that are critical to the modeling and di- 
agnostics of hot, dense plasmas consist almost entirely of 
theoretical calculations performed in LS'-coupling. High- 
resolution absolute photoionization cross-section mea- 
surements are therefore necessary to benchmark the va- 
lidity of this term-resolved data. 

Due to the shortcomings of solar studies, abundance 
determinations shed more light on the argon abundance 
at near-solar metallicity. Such an approach has been sug- 
gested by recent studies on a similar case of neon in B- 
type stars in the Orion nebula [6jX ^^^^ i^ the interstellar 
medium (ISM) toward the Crab Nebulae [S]. Photoion- 
ization cross-section data on this C£-like ion are required 
to accurately determine elemental abundances, for ex- 
ample in the whole stellar sample of B-type stars where 
the A = 442.6 nm and A = 443.0 nm lines of Ar"*" are 
vividly seen in the B main-sequence stars of the Orion 
nebulae j^. The non-LTE abundance calculations were 
based on photoionization cross section data taken from 
the OPACITY PROJECT [3] calculated in LS'-coupling. 
It has been shown for a variety of ions, that the OPAC- 
ITY PROJECT data has severe limitations when com- 
pared to high resolution experimental data |10llllj . 

Photoionization phenomena are bound-free processes 
characterized by interfering channels. In direct photoion- 



ization of a positive ion, the cross section rises from 
zero as a step-function at the ionization threshold en- 
ergy and falls off monotonically with increasing photon 
energy. This process leads immediately to the production 
of Ar^"*" and a free electron 
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Superimposed upon this direct photoionization cross sec- 
tion are series of resonances occurring at discrete photon 
energies corresponding to the excitation of autoionizing 
states in the singly charged Ar"*" ion. Interference be- 
tween the indirect and direct photoionization channels 
produces characteristic Fano-Beutler line profiles for the 
resonances [12', providing further insight into the elec- 
tronic structure of the ion and the dynamics of the pho- 
toionization process. 

A thorough compilation of merged-beams photoion- 
ization studies has been presented in a topical review 
by Kjeldsen |T3] and a previous review of West |14j . 
The technique of photoion-yield spectroscopy using syn- 
chrotron radiation was pioneered by Lyon and collabora- 
tors [15 at the Daresbury synchrotron in the U.K., and 
involves merging ion beams accelerated to keV energies 
with monochromatic, tunable beams of synchrotron ra- 
diation. Photoions produced by their interaction over 
common paths of tens of centimeters may be magneti- 
cally or electrostatically separated from the parent ion 
beam. The directed ion beam and negligible momentum 
transfer permit a complete collection of the photoions. If 
the energy of the ion beam is not too high and ultra-high 
vacuum conditions prevail, the background produced by 
stripping collisions of the primary ion beam with residual 
gas is manageable, and may be subtracted by chopping 
the photon beam. The merged-beams technique is partic- 
ularly amenable to absolute cross-section measurements 
[16]. 

Until recently, because of the relatively low density pre- 
sented by a space-charge-limited ion beam (^10^ cm^'^), 
data for ions were limited mainly to singly-charged al- 
kali, alkaline earth and transition metal ions with large 
photoionization cross sections [14] . The advent of third- 
generation synchrotron radiation sources (3GLS) with 
insertion devices increased the photon beam intensities 
available to experimenters by nearly three orders of mag- 
nitude, making possible the study of photoionization of 
ions with unprecedented sensitivity and spectral resolu- 
tion. This has permitted measurements with ions of the 
lighter and more astrophysically abundant elements, e.g.; 
C+ [17], 0+ [I8l[19], C2+ |2Q] and more recently C^+ [H], 
N'^"'" and O'*"'" [55] , which are helping to refine theoretical 
descriptions of the photoionization process. 

As previously indicated the properties of Argon ions 
are of importance in atmospheric and astrophysical plas- 
mas and are widely used in laboratory plasmas as di- 
agnostic impurities [23J [M|. Properties of Argon ions 
also find many applications in plasma etching and lasers 



|25| . In the present study high-resolution absolute pho- 
toionization cross-section measurements and theoretical 
calculations are presented for the Ar+ ion at photon en- 
ergies ranging from the photoionization threshold to 60 
eV. 



II. EXPERIMENT 

The experiment was conducted on undulator beamline 
10.0.1.2 of the Advanced Light Source (ALS) at Lawrence 
Berkeley National Laboratory. The ion-photon-beam 
(IPB) endstation is based on the merged-beams tech- 
nique and photoion spectroscopy using tunable syn- 
chrotron radiation. The apparatus and methodology 
have been described in detail previously for measure- 
ments of photoionization cross sections for Ne+ ions [TD] . 
Initial measurements of photoionization cross sections for 
0+ ions ^ISj and C^+ ions ^20j using the IPB endsta- 
tion were reported in the literature. The ALS group has 
recently reported detailed experimental and theoretical 
studies on the C^-like Calcium (Ca^+) ion [26, in addi- 
tion to a variety of other singly and multiply charged ions, 
particularly of astrophysically interest, using this same 
IPB endstation; He-like; Li+ [27, 28 , Li-like; B^+ [29J, 
Be-like; C^+ [3U], B+ [31], C^+, N^+ and 0"*+ [211131], 
B-like; C+ [33], F-like; Ne+ [TU], N-like ions; 0+ [18], 
F^"*" and Ne^"*" |TT]. These high resolution experimental 
results enable theoretical studies to be benchmarked, so 
that the resulting data can be used with confidence in 
various astrophysical applications that utilize modeling 
codes such as Cloudy [31135] and XSTAR [35]. 



A. Photon beam 

The photon beam was produced by a 10-cm period 
undulator installed in the 1.9 GeV electron storage ring 
to serve ALS Beamline 10.0.1. A grazing-incidence 
spherical-grating monochromator delivered a highly coUi- 
mated photon beam of spatial width less than 1 mm and 
divergence less than 0.5°. The beamline produces a pho- 
ton fiux of 5 X 10^^ photons per second in a bandwidth 
of 0.01% at an energy of 40 eV and three gratings cover 
the energy range 17-340 eV. Spectral resolving powers 
E/AE as high as 40,000 are available at reduced pho- 
ton fiux. All of the measurements with Ar+ ions were 
carried out with a gold-surfaced spherical grating ruled 
at 380 lines/mm. The photon energy was scanned by 
rotating the grating and translating the exit slit of the 
monochromator while simultaneously adjusting the un- 
dulator gap to maximize the beam intensity. The spec- 
tral resolution was pre-selected by adjusting the entrance 
and exit slits of the monochromator. The photon fiux was 
measured by an absolutely calibrated silicon X-ray pho- 
todiode, and was typically 2-3x10^^ photons/second at 
a nominal spectral resolving power of 2,000. The analog 
output from a precision current meter was directed to 



a volt age-to- frequency converter, which provided a nor- 
malization signal to a personal-computer-based data ac- 
quisition system. The photon beam was time modulated 
(mechanically chopped) at 0.5 Hz using a stepping-motor 
controlled paddle to separate photoions from background 
produced by stripping of the parent ion beam on residual 
gas in the ultra-high vacuum system. The photon energy 
scale was calibrated using measurements ^IS; of the well- 
known 0+ ground-state ('^5') and metastable-state {^P 
and ^D) energy thresholds, allowing for the Doppler shift 
due to the ion motion in the laboratory frame. The abso- 
lute uncertainty in the photon energy scale is estimated 
to be ±10 meV. 



B. Ion beam 

The ^"^Ar^ ions were produced by the Cuernavaca Ion 
Gun Apparatus (GIGA) in a hot-filament, low-pressure 
discharge- type ion source, and accelerated to an energy of 
6 keV. The ion beam was focused by a series of cylindrical 
electrostatic einzel lenses and the Ar+ ion beam selected 
by a 60° analyzing magnet with a mass resolving power of 
100. The ion beam trajectory and its cross-sectional area 
were defined by adjustable horizontal and vertical beam 
slits located downstream of the analyzing magnet. The 
coUimated Ar"*" beam was typically a few mm in diameter, 
with a current in the range 500-650 nA. Three stages of 
differential pumping with turbomolecular pumps and a 
cryopump assured a downstream vacuum in the lO"^*' 
torr range with beams present. 



C. Merger and beam interaction region 

A pair of 90° spherical-sector bending-plates merged 
the ion beam onto the axis of the counter-propagating 
photon beam. Fine-tuning of the overlap of the beams 
was achieved with two sets of mutually perpendicular 
electrostatic steering plates mounted immediately be- 
fore the merger plates. A cylindrical einzel lens focused 
the beam in the center of the interaction region, which 
consisted of an isolated stainless-steel-mesh cylinder to 
which an electric potential (typically +2 kV) was ap- 
plied, thereby energy-labeling photoions produced in this 
region. Series of entrance and exit apertures accurately 
defined the effective length (29.4 cm) of the interaction 
region. Two-dimensional intensity distributions of both 
beams were measured by commercial rotating-wire beam 
profile monitors installed just upstream and downstream 
of the interaction region, and by a translating-slit scanner 
located in the middle of the region. The profile monitors 
permitted the positions and spatial profiles of the two 
beams to be continuously monitored on an oscilloscope 
while tuning the beams. Two 500 1/s mag-ion pumps as- 
sured ultra-high vacuum conditions in this region when 
the photon and ion beams were present. 



D. Demerger and ion charge analyzer 

A 45° dipole analyzing magnet located downstream of 
the interaction region demerged the beams and separated 
the Ar^"*" products from the parent Ar+ beam, which was 
collected in an extended Faraday cup. The magnetic field 
was set such that the Ar^+ product ions passed through 
an aperture in the back of the Faraday cup. A spherical 
90° electrostatic deflector directed them onto a stainless 
steel plate biased at -550V, from which secondary elec- 
trons were accelerated and detected by a microsphere- 
plate electron multiplier used in a pulse-counting mode. 
The deflection planes of the demerger magnet and this 
spherical deflector were orthogonal, permitting the Ar^+ 
products to be swept across the detector in mutually 
perpendicular directions, providing a diagnostic of their 
complete collection. To this end, a cylindrical einzel lens 
located downstream of the interaction region provided 
a further diagnostic, but was found to be unnecessary 
and was turned off during the measurements. The abso- 
lute efficiency of the photoion detector was calibrated in 
situ using an averaging sub-femtoampere meter to record 
the Ar^+ photoion current, which was then compared to 
the measured photoion count rate. The primary Ar+ 
ion beam current was measured by a precision current 
meter, whose analog output was directed to a voltage-to- 
frequency converter, providing a normalization signal to 
the data acquisition system. 



III. THEORY 

For comparison with high-resolution measurements 
such as those at ALS, state-of-the-art theoretical meth- 
ods are required using highly correlated wavefunctions. 
In addition, relativistic effects are required since fine- 
structure effects can be resolved. As metastable states 
are populated in the Ar+ ion beam, additional theoretical 
calculations were required to assist in the determination 
of their fraction. Similar features have previously been 
demonstrated by past detailed experimental and theoret- 
ical photoionization studies on a number of simple and 
complex ions, as outlined in section II. 

Previous photoionization cross-section calculations for 
Ar+ were carried out in L5'-coupling for the OPAGITY 
PROJEGT [3 . To provide a benchmark comparison 
and a guide to the experiment, photoionization cross- 
section calculations were performed in intermediate cou- 
pling (Breit-Pauli) using the R-matrix method [37 . The 
inclusion of relativistic effects is necessary due to the high 
energy resolution of the measurements (nominally 10 
meV) . This theoretical approach optimizes use of limited 
experimental beam time, since coarse energy scans may 
be used in resonance-free regions and fine scans at high 
resolution may be concentrated in energy ranges where 
dense resonance structure is predicted. In the theoretical 
work twenty-four LS terms of the product Ar^+ ion were 
used, arising from the configurations; ls^2s^2p®3s^3p'*. 



Is22s22/3s3/, ls'^2s^2p^is^-ip^i(f and ls^2s'^2p^-ip^ . 
The high resolution spectral scan measurements show 
resonance peaks with an energy separation of magni- 
tude comparable to the fine-structure splitting of the 
ls^2s'^2p^is'^ip^ "^P^l^ 1/2 levels of Ar+. 

In the senii-relativistic R-matrix calculations per- 
formed for this complex, the 24 LS terms give rise to 
48 LSJ levels for the product Ar^+ ion core, all of which 
were included in the close-coupling photoionization cross- 
section calculations. For the Ar^+ product ion the orbital 
basis set used to describe the 24 target LS terms was 
limited to n=3 in constructing the multi-configuration 
interaction wave functions used in the calculations. Pho- 
toionization cross sections were determined for both the 
^P^,2 ground and the metastable ^-Pwj initial states of 
the Ar+ ion due to the population of metastable states 
in the primary ion beam used in the measurements. 

An efficient parallel version |3S] of the R-matrix pro- 
grams |39fHT] was used to determine all the photoioniza- 
tion cross sections. The scattering wave functions were 
generated by allowing for triple-electron promotions out 
of the base configuration 3s^3p^ of Ar+ into the n=3 or- 
bital set employed for the photoionization cross sections 
for both the ground "^P^i^ and the metastable '^P1i2 iiii- 
tial states of Ar+ in intermediate coupling {LSJ) . 

All the scattering calculations were performed with 
twenty continuum basis functions and a boundary radius 
of 7.0 Bohr radii. In the case of the '^P^i2 initial ground 
state, the dipole selection rule requires the dipole transi- 
tion matrices; 3/2° -^ l/2^3/2^5/2^ to be calculated, 
whereas for the metastable '^-Pwj inital state, only the 
dipole matrices for the transitions; 1/2° — > 1/2^^,3/2'^ 
are required. The Hamiltonian matrices for the 1/2°, 
3/2°, 5/2^^, 3/2*^ and 1/2"^ symmetries were calculated 
over the entire range of LS partial waves that contribute 
to these Jtt symmetries. The Ar'^^ theoretical energies 
were shifted to experimental values in the diagonaliza- 
tion procedure of the appropriate Hamiltonians so that 
the energy positions of resonances relative to each ioniza- 
tion threshold would be improved |42l I43| . This energy 
shift was less than 1% for all the energy levels involved. 
For both the ^^^3/2 ground and the metastable ^Pii2 ini- 
tial states, the electron- ion collision problem in the outer 
region was then solved with the parallel R-matrix codes 
using a fine energy mesh of 0.68 /ieV (5 x 10~* Ryd- 
bergs). This ensured that the fine resonance structure in 
the respective photoionization cross sections was fully re- 
solved, particularly in the near-threshold energy region. 
The absolute theoretical photoionization cross sections 
were convoluted with a Gaussian with a FWHM corre- 
sponding to the experimental resolution (nominally 10 
meV) and statistically weighted in order to compare di- 
rectly with the experimental data. 

The multi-channel R-matrix QB technique (applica- 
ble to atomic |44|l45| and molecular complexes |46]) was 
used to determine the resonance parameters. In the QB 
method, the resonance width F is determined from the 



inverse of the energy derivative of the eigenphase sum 5 
at the position of the resonance energy Er via 



F = 2 



dE 



= 2[<5']^L 



E = Er 



(2) 



E=E,^ 



This approach exploits the analytical properties of the R- 
matrix to obtain the energy derivative of the reactance 
(K) matrix and has previously proved very successful in 
locating fine resonance features in the cross sections for 
multi-channel electron scattering. 



IV. RESULTS 

A. Absolute cross-section measurements 

Absolute measurements of photoionization cross sec- 
tions were performed at a number of discrete photon 
energies where there are no resonant features in the 
photoion-yield spectrum. At each such photon energy 
(hiy), the value of the total absolute photoionization cross 
section crpi in cm^ was determined from experimentally- 
measured parameters: 
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(3) 



where R is the photoion count rate [s"^], q is the charge 
state of the parent ion, e = 1.60 x 10~^^ C, Vi is the ion 
beam velocity [cm/s] , e is the responsivity of the photo- 
diode [electrons/photon] , /+ is the ion beam current [A] , 
/ "^ is the photodiode current [A] , il is the photoion collec- 
tion efficiency, 61 is the pulse transmission fraction of the 
photoion detection electronics (determined by the pulse- 
discriminator setting), A is the measured absolute pho- 
toion detection efficiency, and the beam overlap integral 
jF{z)dz defines the spatial overlap of the photon and 
ion beams along the common interaction path in units 
of cm~^. The propagation direction of the ion beam is 
defined as the z-axis. At each of the three positions z^ at 
which beam intensity profiles were measured, the form 
factor F{zi) was determined by the following relation. 



F{Z^) 
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(4) 



The absolute photoion detection efficiency was deter- 
mined in situ by comparing the count rate of Ar^+ ions on 
the detector with the current measured using a calibrated 
averaging sub-femtoammeter. In Table IT] experimental 
parameters are given for a typical absolute cross-section 
measurement at a photon energy of 32.5 eV. Such mea- 
surements were repeated several times at each photon 
energy. 

Table [IT] presents the relative, absolute and the total 
uncertainties (their quadrature sums) in the photoioniza- 



TABLE I: Values of experimental parameters for a typical 
Ar+ absolute photoionization cross-section measurement at 
32.5 eV. 



Parameter 



Value(s) 



Ion Beam Energy 

Ion Beam Current, J+ 

Photon Energy, hu 

Photodiode current, P 

Photon flux 

Interaction bias voltage, Vi„t 

Ion interaction velocity, Vi 

Ar^"*" signal rate, R 

Ar'^"'" background rate 

Form Factors: F(zi), F{z2), F{zz] 

Photodiode responsivity, e 

Merge-path length, L 

Pulse transmission fraction, 5i 

Photoion collection efficiency, Q. 

Photoion detection efficiency, A 

Cross Section, a 



6.0 keV 

620.0 nA 

32.4 eV 

19.7 M 

5.6x10" photons/s 

+2.0 kV 

1.39x10^" cm/s 

7012.5 s'^ 

45.2 s-^ 

41.22, 28.22, 19.17 cm" 

7.81 electrons/photon 

29.4 cm 

0.75 

1.00 

0.216 

10.95 Mb 

(1Mb = 10'" cm^) 



TABLE II: Systematic uncertainties in absolute cross-section 
measurements estimated at 90% confidence level. 



Source 



Relative Absolute Total 



Counting Statistics 
Photoion Detector Efficiency 
Photoion Collection Efficiency 
Pulse Counting Efficiency 
Primary Ion Collection Efficiency 
Ion Current Measurement 
Photodiode Responsivity 
Photodiode Current Measurement 
Beam Profile Measurement 
Beam Overlap Integral 
Interaction Length 



2% 


- 


2% 


- 


10% 


10% 


2% 


2% 


3% 


- 


3% 


3% 


- 


2% 


2% 


- 


2% 


2% 


5% 


10% 


11% 


2% 


2% 


3% 


3% 


7% 


8% 


10% 


7% 


12% 


- 


2% 


2% 



Quadrature Sum 



12% 



18% 22% 



tion cross-section measurements estimated at 90% confi- 
dence level (two standard deviations on statistical uncer- 
tainties) or at an equivalent confidence level for system- 
atic uncertainties. Table |III| presents the absolute pho- 
toionization cross-section values that were determined at 
the selected photon energies. The data from Table [III| are 
plotted in Fig.jTjwith error bars representing the total un- 
certainties. Photon energy scans taken over the energy 
range 27.4 - 60 eV at a resolution of 10 meV with 1.5 
meV steps were normalized to these absolute measure- 
ments by fitting a polynomial to their measured ratio as 
a function of photon energy. These data are presented in 
Fig. [T] as small triangles. 

Fig. |2] presents the ALS measured photoionization 
cross section in the region of the ^P3/2 ground-state and 
^Pi/2 metastable-state ionization threshold energies at 
27.630 eV and 27.452 eV, respectively according to those 
given by the NIST tabulations |32]- The ratio of the mag- 
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FIG. 1: (Color online) Overview of experimental cross- 
section measurements and theory for Ar^ (shaded area) at 
a resolution of 10 meV over an energy range from the ioniza- 
tion threshold to 60 eV; (a) cross sections on a logarithmic 
scale, (b) cross section on a linear scale. The absolute cross- 
section measurements are indicated by open triangles with 
error bars. The curve represents a 48-level close-coupling 
Breit-Pauli R-matrix calculation performed in intermediate 
coupling. For comparison, the theory curve is a statistically 
weighted sum of cross sections for photoionization from the 
ground and metastable states and has been convoluted with 
a FWHM Gaussian of 10 meV. 



TABLE III: Measured values of the total absolute cross sec- 
tions for photoionization of Ar^ . 



Energy (eV) 
(eV) 



Cross Section 
(10-" cm^) 



Total Uncertainty 
(10-"cm2) 



29.5 
32.5 
36.0 
42.0 
47.0 
49.0 
56.0 
59.0 



21.4 
11.6 
6.0 
1.2 
0.9 
1.0 
1.3 
1.4 



4.6 
2.5 
1.3 
0.3 
0.2 
0.2 
0.3 
0.3 



nitudes of the two threshold cross-section steps is 2.06 
± 0.10, consistent with a statistical population (propor- 
tional to 2 J -I- 1) of ground- and metastable-state ions 
in the primary ion beam. The present measurements 
are therefore considered to correspond to a sum of the 
^^3/2 ground-state cross section multiplied by 2/3 and 
the metastable-state ^Pi/2 cross section multiplied by 
1/3. The presence of the metastable component adds 
complexity to the observed resonance structure and to its 
interpretation. A very small offset of the measured cross 
section from zero below the ^Pi /2 threshold is attributed 
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FIG. 2: (Color online) High resolution PI cross-section mea- 
surements from the ALS at 10 meV energy resolution in the 
photon energy region of the ^^3/2 ground-state and ■^Pi/2 
metastable-state ionization thresholds of Ar^ at 27.630 eV 
and 27.452 eV, respectively. The small, non-zero cross section 
at energies below both thresholds is attributed to higher-order 
radiation in the photon beam. Theoretical results are from an 
intermediate coupling Breit-Pauli R-matrix calculation, sta- 
tistically weighted for the ground and metastable state. The 
theoretical results have been convoluted with a Gaussian of 
10 meV FWHM to simulate the photon energy bandwidth of 
the experiment. 
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FIG. 3: (Color online) Experimental results from the ALS 
for the PI cross sections of Ar^ for photon energies ranging 
from 27-32 eV with a photon energy resolution of 10 meV. 
In the figure, the members of Rydberg series of resonances 
originating from the ^P^i2 metastable (|) and ^P^i2 ground 
states (A) of Ar^ are identified. The resonance assignments 
were determined from the NIST tabulated levels for Ar II and 
Ar III |42]. 



to the presence of higher-order radiation from the undu- 
lator, which is dispersed by the grating and estimated to 
comprise 6% of the total photon flux in this energy region. 
In this threshold region, a detailed comparison of the 
ALS experimental results with theoretical results (from 
a Breit-Pauli calculation performed in intermediate cou- 
pling) show excellent agreement. In order to make a di- 
rect comparison with experiment in this energy region the 
theoretical results were convoluted with a FWHM gaus- 
sian of 10 meV and statistically weighted over the ground 
and metastable states. The two strong resonance fea- 
tures observed in the threshold energy region of the ALS 
photoionization experimental measurements, (located at 
27.481 eV and 27.658 eV respectively; see Fig [2] ) re- 
sulting in photoexcitation of the ?>s^?>p^{}D2)Qd Auger 
state, (from either the ground state or metastable states 
of Ar+), are reproduced suitably well by the present theo- 
retical Breit-Pauli calculations performed in intermediate 
coupling. 

The peak heights of the stronger calculated resonances 
exceed the measured values, and the experimental line- 
shapes of the resonances at 27.48 eV and 27.66 eV are 
non-Gaussian. This suggests a possible saturation effect 
due to the dead time of the photoion counting system at 
the higher signal counting rates. The measured ratios of 
the strengths of these two resonance further suggests that 
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FIG. 4: (Color online) Experimental results from the ALS for 
the PI cross sections of Ar^ for photon energies ranging from 
27-32 eV with a photon energy resolution of 10 meV compared 
with theoretical results from a 48-state Breit-Pauli R-matrix 
calculation. The theoretical results have been convoluted with 
a Gaussian of FWHW of 10 meV and statistically averaged 
over the ground and metastable states. 
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FIG. 5: (Color online) Eigenphase sum S and its derivative 
S' for each 2J^ contributing to the total PI cross section a 
originating from the initial Ar"*'(^P°/2) metastable state. The 
PI cross section has been convoluted with a FWHM Gaussian 
of 10 meV to simulate the photon experimental bandwidth. 



FIG. 7: (Color online) Linewidths (meV) for the domi- 
nant 3s^3p'^{^So)ns,nd resonances series converging to the 
^So threshold of the product Ar^"*" ion originating from the 
initial Ar"''(^Pj°/2) metastable state. 
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the ground-state fraction in the primary ion beam may 
have been closer to 0.73 than the value of 0.67 based on 
statistical -weights that "was used to scale the theoretical 
results for comparison "with experiment. 

In Fig.[3]"we present the experimental data taken at 10 
meV in the photon energy region from threshold to ap- 
proximately 32 eV. As can been seen from Fig. [3] there is 
extremely rich resonance structure evident in photoion- 
ization cross section data. The assignments of the var- 
ious Rydberg resonance series are given based on the 
NIST tabulated levels for Ar II and Ar III 02]. Only the 
prominent resonance series are analyzed in detail. Fig. HI 
illustrates the same energy range but includes the com- 
parison "with the theoretical results obtained from the in- 
termediate coupling Breit-Pauli calculations. Here again 
theory has been convoluted "with a FWHM Gaussian of 
10 meV and statistically "weighted over the ground and 
metastable states to make a direct comparison "with ex- 
periment in this energy region. Once again it is seen 
there is remarkable agreement bet"ween theory and ex- 
periment, both on the absolute cross section scale and 
on the photon energy scale. 



FIG. 6: (Color online) Linewidths (meV) for the dominant 
3s^3p'^{^D2)ns,nd resonances series converging to the ^ D2 
threshold of the product Ar^""" ion originating from the initial 
At+CP^^2) metastable state. The 3s^3p''(^D2)ng Rydbergy 
resonance series have linewidths very much less than 10~^ 
meV and have not been analyzed. 



B. Resonances 

In the photon energy range 27 eV to 32 eV a "wealth of 
resonance structure is observed in the PI cross sections as 
illustrated in Fig. 3. We note that belo"w the ^D2 thresh- 
old of the product AP+ ion, 3s^3p^{^D2)ns "^D, 
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FIG. 8: (Color online) Eigenphase sum S and its derivative 
S' for each 2J^ contributing to the total PI cross section a 
originating from the initial Ar^(^P3/2) ground state. The PI 
cross section has been convoluted with a FWHM Gaussian of 
10 meV to simulate the photon experimental bandwidth. 



FIG. 10: (Color online) Linewidths (meV) for the domi- 
nant 3s^3p*(^5'o)ns,nd resonances series converging to the 
^So threshold of the product Ar^""" ion originating from pho- 
toionization of the initial Ar^(^P3/2) ground state. 
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FIG. 9: (Color online) Linewidths (meV) for the domi- 
nant 3s^3p'^{^ D2)ns,nd resonances series converging to the 
^D2 threshold of the product Ar^^ ion originating from pho- 
toionization of the initial Ar^(^P3/2) ground state. The 
3s^3p'^{^D2)ng Rydbergy resonance series have linewidths 
very much less than 10~^ meV and have not been analyzed. 



Ss 3p { D2)nd Si/2 dominant Rydberg resonance series 
may occur from the angular momentum coupling rules. 

At higher photon energies 3s^3p'^{^So)nd ^-03/2.5/2 £^nd 
3s^3p^("'^S'o)ns ^5*1/2 Rydberg resonance series converg- 
ing to the ^Sq threshold of the corresponding product 
Ar^+ ion are observed. The resonance series found are 
illustrated in Figs. 5-10. For resonances converging to 
both the ^D2 and ^Sq thresholds of the Ar^+ ion only 
the dominant resonances series are analyzed. 

To assist in the identification of the dominant reso- 
nance features observed in the photoionization spectra a 
theoretical analysis was carried out based on the widely 
used QB technique for determining resonance parameters 
in atomic and molecular systems |iil - liS] . This technique 
exploits the properties of the R-matrix in multi-channel 
scattering and is particularly helpful for locating and de- 
termining the properties of narrow resonances found in 
the spectra. 

The relationship between the principal quantum num- 
ber n, the effective quantum number v and the quantum 
defect ^ for an ion of effective charge Z is given by j/ = 
n - /z where the resonance position e^ can be determined 
from Rydberg 's formula 



zv-^ 



(5) 



where Coo is the resonance series limit. 

Prominent Rydberg resonance series observed in the 
experimental photoionization spectra have been analyzed 
and assigned according to their quantum defects. The 
NIST tabulations were used as a helpful guide in the as- 
signment. The principal quantum number n. resonance 



9 



energy in eV, quantum defect /x and the resonance au- 
toionizing linewidth F (nieV) of the first few members of 
the dominant Rydberg series originating from the ^P°/2 
metastable and the '^P'^/2 ground state of the Ar+ ion 
are listed in the Tables IV, V and VI. The results for 
the appropriate resonance parameters for the dominant 
Rydberg resonance series found in the spectra of the 
PI cross section for the ground and metastable states 
of Ar+. lying below the ^D2 and ^Sq thresholds of the 
product Ar^+ ion are listed in Tables IV, V and VI. As 
may be seen from the results of the QB analysis, many 
members of the various Rydberg series have resonance 
linewidths smaller than the experimental photon band- 
width of 10 meV, therefore making their detection ex- 
tremely difficult. Finally we note that interloping res- 
onances 3s^3p^(^S'o)6s, 5d are present in the spectrum 
below the 3s'^3p'^{^D2) threshold (as clearly illustrated 
in figures 6 and 9) disrupting the regular pattern of the 
Rydberg series. 



V. SUMMARY 

Absolute photoionization cross sections have been 
measured for a statistical mixture of ground-state and 
metastable Ar+ in the photon energy range w 27 eV - 
60 eV with a nominal energy resolution of 10 meV. The 
cross section in this energy region is characterized by a 
wealth of resonance features due to indirect photoion- 
ization. The dominant features were spectroscopically 
assigned to Rydberg series and quantum defects were de- 
termined. A conservative estimate for the accuracy of the 
resonance energies is ± 10 meV. The natural linewidths 
were determined from the R-matrix method which we 
would estimate to be accurate to 10 %. Over the en- 
tire energy range investigated, the theoretical results for 
the photoionization cross sections on this C^-like system 
obtained from an intermediate coupling Breit-Pauli R- 
matrix approach are in excellent agreement with the high 



resolution measurements. The main thrust of this joint 
investigation is the coupling of high experimental sensi- 
tivity and photon energy resolving power with state-of- 
the-art theoretical predictions. The excellent agreement 
between theory and experiment over the entire energy 
range gives added confidence in both results. The pho- 
toionization cross sections from the present study on this 
C£-likc ion are suitable for inclusion in state-of-the-art 
photoionization modeling codes such as Cloudy |34l [55] 
and XSTAR [36J , which are used to numerically simulate 
the thermal and ionization structure of ionized astrophys- 
ical nebulae. 



VI. ACKNOWLEDGMENTS 

The experimental work was supported by the Office of 
Basic Energy Sciences, Chemical Sciences, Geosciences 
and Energy Biosciences Division, of the U.S. Department 
of Energy under grants DE-FG03-00ER14787 and DE- 
FG02-03ER15424 with the University of Nevada, Reno; 
by the Nevada DOE/EPSCoR Program in Chemical 
Physics and by CONACyT through the CCF-UNAM, 
Cuernavaca, Mexico. A. A. and M. M. S'A. acknowledge 
support from DGAPA-UNAM-IN 113010 (Mexico) and 
CNPq (Brazil), respectively. C. P. B. was supported by 
U.S. Department of Energy grants through Auburn Uni- 
versity. B. M. McL. thanks the Institute for Theoretical 
Atomic and Molecular Physics (ITAMP) for their hospi- 
tality and support under the visitors program. ITAMP 
is supported by a grant from the U.S. National Science 
Foundation to Harvard University and the Smithsonian 
Astrophysical Observatory. The computational work was 
carried out at the National Energy Research Scientific 
Computing Center in Oakland, CA and on the Tera-grid 
at the National Institute for Computational Sciences in 
Knoxville, Tennessee, which is supported in part by the 
U.S. National Science Foundation. 



[1] J. N. Bregman and J. P. Harrington, Astrophys. J. 309, 

833 (1986). 
[2] I. Hofmann, Laser Part. Beams 8, 527 (1990). 
[3] W. Cunto, C. Mendoza, F. Ochsenbein and C. J. Zeippen, 

Astron. Astrophys. 275, L5 (1993). 
[4] C. A. Iglesias and F. J. Rogers, Astrophys. J. 464, 943 

(1996). 
[5] M. A. Bautista and T. R. Kallnian, Astrophys. J. Suppl. 

Ser. 134, 139 (2001). 
[6] K. Cunha, I. Hubeny and T. Lanz, Astrophys. J. 647, 

L143 (2006). 
[7] C. Esteban, M. Peimbert, J. Garci'a-Rojas, M. Ruiz, A. 

Peimbert and M. Rodriguez, Mon. Not. Roy. Astro. Soc. 

355, 229 (2004). 
[8] J. S. Kaastra, C. P. de Vries, E. Constantini and J. W. 

A. den Herder, Prog. Theo. Phys. 169, 204 (2007). 



[9] T. Lanz, K. Cunha, J. Holtzman and 1. Hubeny, Astro- 
phys. J. 678, 1342 (2008). 

[10] A. M. Covington, A. Aguilar, L R. Covington, M. F. 
Charaibeh, G. Hinojosa, C. A. Shirley, R. A. Phaneuf, 
L Alvarez, C. Cisneros, I. Dominguez-Lopez, M. M. 
Sant'Anna, A. S. Schlachter, B. M. McLaughlin and A. 
Dalgarno, Phys. Rev. A 66, 062710 (2002). 

[11] A. Aguilar, E. D. Emmons, M. F. Charaibeh, A. M. Cov- 
ington, J. D. Bozek, G. Ackerman, S. Canton, B. Rude, 
A. S. Schlachter, C. Hinojosa, I. Alvarez, C. Cisneros, B. 
M. McLaughhn and R. A. Phaneuf, J. Phys. B: At. Mol. 
& Opt. Phys. 38, 343 (2005). 

[12] U. Fano and J. W. Cooper, Phys. Rev. A 137, 1364 
(1965). 

[13] H. Kjeldsen, J. Phys. B: At. Mol. & Opt. Phys. 39, R325 
(2006). 



10 



TABLE IV: Principal quantum numbers n, resonance energies (eV), quantum defects /i and autoionization linewidths [r(meV)], 
determined by the QB method for the dominant 3s'^3p^{^So)ns,nd series observed in the photoionization spectra originating 
from the 3s^3p^ '^Pi/2 metastable state of Ar""". 







Rydberg Series 




Width (r) 


Rydberg Series 




Width (r) 






3s^3p-^{^D2)ns,nd 




3s^3p'^CD2)ns,nd 


3s 3p C' So)ns,nd 




3s 3p ( So)ns,nd 


Initial State 


n 


Energy (eV) 


^J■ 


Energy(meV) 


Energy (eV) 


M 


Energy (meV) 








J=l/2 










Ss^Sp^ 2p°/2 


|7s| 


- 


- 


- 


29.645 


1.692 


1.36 




8 


27.830 


1.673 


103.4 


30.209 


1.691 


0.39 




9 


28.173 


1.682 


64.8 


30.558 


1.690 


0.14 




10 


28.401 


1.688 


43.0 


30.789 


1.689 


0.06 




11 


28.561 


1.692 


30.0 


30.949 


1.688 


0.03 




12 


28.677 


1.695 


21.8 


31.065 


1.688 


0.01 




13 


28.763 


1.698 


16.3 


31.151 


1.688 


< 10"^ 




14 


28.829 


1.699 


12.5 


31.217 


1.688 


<io-2 




15 


28.881 


1.701 


9.8 


31.269 


1.688 


< 10"^ 




16 


28.923 


1.702 


7.8 


31.311 


1.687 


<io-2 



29.189 



31.596 



3s23p^ 2p° 



/2 



J=l/2 



[6dl 


- 


- 


- 


7 


27.830 


0.521 


0.16 


8 


28.211 


0.540 


0.11 


9 


28.426 


0.554 


0.07 


10 


28.576 


0.574 


0.05 


11 


28.691 


0.546 


0.04 


12 


28.773 


0.562 


0.03 


13 


28.837 


0.568 


0.02 


14 


28.887 


0.572 


0.02 


15 


28.928 


0.574 


0.02 



29.189 



J = 3/2 



3s23p5 2p°/2 


[7s] 


- 


- 






8 


27.821 


1.692 


0.44 




9 


28.171 


1.691 


0.22 




10 


28.401 


1.690 


0.13 




11 


28.561 


1.689 


0.08 




12 


28.677 


1.688 


0.08 




13 


28.764 


1.688 


0.04 




14 


28.830 


1.688 


0.02 




15 


28.882 


1.687 


0.02 




16 


28.923 


1.687 


0.02 



29.189 



38^3/ 2p° 



J = 3/2 



/2 



[6d] 


- 


- 




7 


27.855 


0.613 


83.4 


8 


28.163 


0.716 


25.7 


9 


28.395 


0.720 


18.1 


10 


(28.590) 


(0.422) 


33.8 


11 


28.672 


0.738 


3.8 


12 


28.760 


0.745 


5.3 


13 


28.827 


0.744 


4.4 


14 


28.879 


0.744 


3.6 


15 


28.921 


0.744 


2.9 



29.615 


0.732 


28.9 


30.185 


0.746 


16.4 


30.540 


0.754 


10.1 


30.775 


0.759 


6.6 


30.938 


0.763 


4.6 


31.057 


0.766 


3.3 


31.145 


0.767 


2.5 


31.213 


0.769 


1.9 


31.266 


0.770 


1.5 


31.308 


0.771 


1.2 



oo" 29.189 

"Series limit from NIST Standard Reference Database 1421. 



31.596 
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TABLE V: Principal quantum numbers n, resonance energies (eV), quantum defects fj, and autoionization linewidths [r(meV)], 
determined from the QB method for the dominant Ss^3p'^{^ D2)ns,nd and 3s^3p*(^S'o)ns,nd series resulting from photoion- 
ization of the 3s^3p^ ^^3/2 ground state of Ar^. 







Rydberg Series 




Width (r) 


Rydberg Series 




Width (r) 






3s^3p'^(}D2)ns,nd 




3s^3p-^{^D2)ns,nd 


3s'^3p^{^So)ns,nd 




3s 3p C'So)ns,nd 


Initial State 


n 


Energy (eV) 


^ 


Energy (meV) 


Energy (eV) 


M 


Energy (meV) 








J = l/2 










3s'3p' 2P3°/2 


|7s| 


- 


- 


- 


29.822 


1.692 


1.36 




8 


28.007 


1.673 


103.4 


30.387 


1.691 


0.39 




9 


28.351 


1.682 


64.8 


30.736 


1.690 


0.14 




10 


28.579 


1.688 


43.0 


30.966 


1.689 


0.06 




11 


28.739 


1.692 


30.0 


31.126 


1.689 


0.03 




12 


28.854 


1.696 


21.8 


31.242 


1.688 


0.01 




13 


28.941 


1.698 


16.3 


31.329 


1.688 


< 10"^ 




14 


29.007 


1.700 


12.5 


31.395 


1.688 


<lQ-2 




15 


29.059 


1.701 


9.8 


31.447 


1.688 


< 10"^ 




16 


29.100 


1.702 


7.8 


31.488 


1.688 


<lQ-2 



29.367 



31.774 



J = l/2 



3s'3p5 'P3°/2 


[6d] 


- 


- 


- 




7 


28.070 


0.521 


0.16 




8 


28.389 


0.540 


0.11 




9 


28.603 


0.554 


0.07 




10 


28.754 


0.574 


0.05 




11 


28.869 


0.546 


0.04 




12 


28.951 


0.562 


0.03 




13 


29.014 


0.568 


0.02 




14 


29.065 


0.572 


0.02 




15 


29.105 


0.574 


0.02 



29.367 



J=3/2 



Ss'Sp" 2p°/2 


[7s] 


- 


- 


- 




8 


27.999 


1.692 


0.44 




9 


28.348 


1.691 


0.26 




10 


28.579 


1.690 


0.13 




11 


28.739 


1.689 


0.08 




12 


28.855 


1.688 


0.06 




13 


28.941 


1.688 


0.04 




14 


29.098 


1.688 


0.03 




15 


29.060 


1.687 


0.02 




16 


29.101 


1.687 


0.02 



29.367 



J=3/2 



Ss'3p''PS/2 


[6d] 


- 


- 


- 




7 


28.033 


0.613 


83.4 




8 


28.341 


0.716 


25.7 




9 


28.573 


0.720 


18.1 




10 


(28.774) 


(0.422) 


33.8 




11 


28.850 


0.738 


3.8 




12 


28.937 


0.745 


5.3 




13 


29.004 


0.744 


4.4 




14 


29.047 


0.744 


3.6 




15 


29.099 


0.744 


2.9 



29.793 


0.732 


28.9 


30.363 


0.746 


16.4 


30.718 


0.754 


10.1 


30.953 


0.759 


6.6 


31.116 


0.763 


4.6 


31.234 


0.766 


3.3 


31.322 


0.768 


2.5 


31.390 


0.769 


1.9 


31.443 


0.770 


1.5 


31.485 


0.771 


1.2 



00" 29.367 

"Series limit from NIST Standard Reference Database 1421. 



31.774 
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TABLE VI: Principal quantum numbers n, resonance energies (eV), quantum defects /j, and autoionization linewidths [F (meV)], 
determined from the QB method for the dominant 3s'^3p^{^D2)ns, nd and 3s^3p''(^S'o)nd series resulting from photoionization 
of the 3s^3p^ '^P-i/2 ground state of Ar+. 







Rydberg Series 




Width (r) 


Rydberg Series 




Width (r) 






3s23p*(^D2)ns,nd 




?,s'^ip'^(^D2)ns,nd 


is^2,p^{^So)nd 




3s23p'*(^S'o)nd 


Initial State 


n 


Energy (eV) 


M 


Energy (meV) 


Energy (eV) 


M 


Energy (meV) 








J=5/2 










Zs^ip^ 2p°/2 


[7s] 


- 


- 


- 










8 


27.999 


1.692 


0.27 










9 


28.348 


1.691 


0.12 










10 


28.579 


1.690 


0.07 










11 


28.739 


1.689 


0.04 










12 


28.855 


1.688 


0.03 










13 


28.941 


1.688 


0.02 










14 


29.008 


1.688 


0.01 










15 


29.060 


1.687 


< 10"^ 










16 


29.101 


1.687 


< 10-2 









29.367 



J=5/2 



3s23p5 '^PS/2 


[6d] 


- 


- 


- 


29.790 


0.735 


25.2 




7 


27.983 


0.729 


2.25 


30.361 


0.749 


14.3 




8 


28.335 


0.737 


1.40 


30.717 


0.757 


8.8 




9 


28.569 


0.742 


0.94 


30.952 


0.763 


5.8 




10 


28.731 


0.745 


0.66 


31.116 


0.766 


4.0 




11 


28.849 


0.748 


0.40 


31.234 


0.769 


2.9 




12 


28.937 


0.750 


0.33 


31.322 


0.771 


2.2 




13 


29.004 
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0.26 


31.390 


0.773 


1.7 




14 


29.057 


0.752 


0.20 


31.443 


0.774 


1.3 




15 


29.099 


0.753 


0.16 


31.485 


0.775 


1.0 
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- 
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